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The BFW model introduced by Bohman, Frieze, and Wormald [Random Struct. Algorithms, 25, 
432 (2004)] and recently investigated in the framework of discontinuous percolation by Chen and 
D'Souza [Phys. Rev. Lett., 106, 115701 (2011)], is studied on the square and simple-cubic lattices. 
In two and three dimensions, we find numerical evidence for a strongly discontinuous transition. In 
two dimensions, the clusters at the threshold are compact with a fractal surface of fractal dimension 
df = 1.49 ± 0.02. On the simple-cubic lattice, distinct jumps in the size of the largest cluster are 
observed. We proceed to analyze the tree-like version of the model, where only merging bonds 
are sampled, for dimension two to seven. The transition is again discontinuous in any considered 
dimension. Finally, the dependence of the cluster-size distribution at the threshold on the spatial 
dimension is also investigated. 

PACS numbers: 64.60.ah, 64.60.al, 89.75.Da 



I. INTRODUCTION 

Percolation is a classical model in Statistical Physics 
which, despite consisting of simple geometrical rules, has 
found application in a broad range of problems ^ [2]. 
In classical (random) percolation, a fraction p of sites or 
bonds are occupied uniformly at random. At a critical 
occupation fraction p c the system undergoes a continuous 
transition from a non-percolating to a percolating state. 
A recent work by Achlioptas, D'Souza, and Spencer [3] 
hinted at the possibility of a discontinuous percolation 
transition by slightly modifying the bond selection rules 
of evolving graphs, coining the term explosive percolation. 
This surprising result led to intensive research efforts 
studying the model on different topologies and dimen- 
sions, yielding ambiguous results regarding the nature of 
the transition [4 18] . More recently, the transition in the 
original model was formally demonstrated to be contin- 
uous in the mean-field limit |19| and there is evidence 
that the same applies to other models [2"Ull23| . Alter- 
native models have been devised and analyzed in detail 
[2~4"rl53| . several of them with clear signs supporting the 
hypothesis of a discontinuous transition. 

In particular, Chen and D'Souza considered the perco- 
lation properties on the random graph of the BFW model 
[341 135) . originally introduced by Bohman, Frieze, and 
Wormald [35] , reporting a discontinuous transition with 
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several giant components in the thermodynamic limit. In 
general, the behavior of a system undergoing a transition 
depends on the dimension [37]. In this paper, we analyze 
the original BFW model on the 2D square lattice and on 
the 3D cubic lattice, as well as the tree-like version up 
to dimension seven. 

The manuscript is organized as follows. In Sec. [IT] we 
discuss the BFW model, adapted to the square lattice, 
along with some definitions. Section [HI| contains the ex- 
tension to three dimensions. The tree-like version is an- 



alyzed in Sec. IV We conclude with the final remarks in 
Sec.E] 



II. BFW MODEL ON THE SQUARE LATTICE 

The BFW model has recently been analyzed by Chen 
and D'Souza [34] [35] on the random graph, by extend- 
ing the algorithm initially introduced by Bohman, Frieze, 
and Wormald [35]. Here, we consider this model on a 
square lattice of linear size L with periodic boundary 
conditions and number of sites N, where N = L 2 . Ini- 
tially, all sites are occupied and all bonds are unoccupied 
(empty), such that there are N clusters of unitary size. 
As in Ref. (3H |3S] , one occupies the bonds according to 
the following procedure. Let u be the total number of 
selected bonds, t the number of occupied bonds, and k 
the stage of the process, initially set to k = 2. The first 
bond is occupied at random, such that t = u = 1. Then, 
at each step u, 

1. One bond, among the unoccupied ones, is selected 
uniformly at random. 

2. I is measured as the maximum cluster size if the 
selected bond is occupied. 
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FIG. 1. (Color online) Fraction of sites in the largest cluster 
Poo (solid red line), stage per site k/N (o), occupied bonds 
per step t/u (dashed upper blue line), and g(k) (dotted lower 
black line) as a function of the fraction of occupied bonds 
p, for the BFW model on the square lattice. As reported in 
Ref. [341 136] for the random graph, k is essentially equal to 
s m „. Therefore the corresponding curves in the figure, Poo 
and k/N, almost overlap. The upper inset shows a close-up 
view of the upper-left corner of the main plot, as indicated by 
the box. The system size is N = 512 2 , the measured values 
are averages over more than 10 5 samples. In the lower inset 
we see the same variables as the main plot, measured for a 
single realization. 



3. If I < k, go to (4). 

Else if t/u > g{k) = 1/2 + y/l/(2k), go to (5). 
Else increment k by one and go to (3). 

4. Occupy the selected bond. 

5. Increment u by one. 

Note that the above procedure, in (1), only samples un- 
occupied bonds, while in Ref. [53H5B] edges are chosen 
uniformly at random from all edges, whether occupied 
or not. In adapting random graph models to the lat- 
tice, we solely sample unoccupied bonds since in classical 
(random) percolation the control parameter is the frac- 
tion of occupied bonds [H 03 [7J [H [2T] . In contrast, for the 
BFW model on the random graph, by sampling all bonds, 
since sufficiently many new links can be accepted asymp- 
totically, more than one stable macroscopic cluster is ob- 
tained [SH|3S]. The procedure is applied iteratively while 
the system evolves from t = 0, corresponding to a bond 
occupation fraction of p = t/ (2N) — 0, to t — 2N where 
all bonds in the system are occupied, i.e., p=l. By con- 
struction, at each stage, the size of the largest cluster can 
be at most s max < k (see Fig. [I]) [36]. To illustrate the 
evolution of this process, we monitor in Fig.[T] the usual 
order parameter for the percolation transition, defined as 
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FIG. 2. (Color online) Fraction of sites in the largest cluster 
Poo as a function of the fraction of occupied bonds p on the 
square lattice for different lattice sizes L (from left to right: 
32, 64, 128, 256, 512, 1024, 2048, 4096, and 8192). 



the fraction of sites in the largest cluster P x = s max /N. 
One observes that the transition is delayed, compared 
to classical percolation where p c = 1/2, and the growth 
of the order parameter at the threshold appears to be 
more pronounced than in the classical case (numerical 
evidence discussed below). In addition, Fig. [I] shows the 
behavior of the stage per site k/N, which grows as the 
system evolves and which imposes a bound on the order 
parameter, P^, < k/N. The fraction of accepted bonds 
t/u and the value of g{k) are also shown as functions 
of the fraction of occupied bonds p. These curves illus- 
trate the mechanism of the BFW model: if the fraction 
of accepted bonds were to drop below g(k), the stage k 
is incremented which decreases g(k) (see insets of Fig. [I]). 
Starting from full acceptance, t/u decreases, approaching 
the asymptotic value of 1/2. The mechanism promotes 
the homogenization of cluster sizes and delays the tran- 
sition. Once percolation occurs, it is more abrupt than 
in the case of uniformly random bond occupation. When 
k/N = = 1, with p < 1, the entire set of unoccupied 
bonds consists of redundant bonds, that is unoccupied 
bonds internal to existing clusters, which can be occu- 
pied without affecting the size of clusters. In Fig. [2] we 
see the fraction of sites in the largest cluster P^ as a 
function of the fraction of occupied bonds p for different 
lattice sizes L. Typical snapshots of the evolving sys- 
tem are shown in Fig.|3j for fractions of occupied bonds 
p = 0.5, 0.9, and 0.95. 

To quantify this observation, we make the following 
considerations. If the largest change J in the order pa- 
rameter Pqo is finite in the thermodynamic limit, a per- 
colation transition is called discontinuous, as opposed to 
continuous transitions where the jump J vanishes for 
L — > oo [38] . Additional ways to distinguish between 
continuous and discontinuous transitions are to consider 
the behavior of the standard deviation of the order pa- 
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FIG. 3. (Color online) Snapshots of the system evolving ac- 
cording to the BFW model, obtained on a square lattice with 
512 2 sites. The following fractions of occupied bonds are 
shown: (a) p = 0.5, (b) 0.9, and (c) 0.95. We note that 
in the t herm odynamic limit p c seems to be equal to unity 
(see Fig. [To]!. 
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where M 2 = J2i s i /N is the second moment of the 
cluster-size distribution and Si is the size of cluster i. Fig- 
ure [4] shows the size independence of the largest change 
(J)j and the order parameter {Poo)] immediately before 
the biggest jump, where the average (-)j means averaging 
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FIG. 4. Order parameter P^ (*) immediately before the 
jump, jump size J (□), and maximum of the second mo- 
ment of the cluster-size distribution M2 (o) as a function 
of the linear system size L. All three quantities are asymp- 
totically independent of the system size. For large system 
sizes, we estimate the following values: J = 0.415 ± 0.005, 
Poo = 0.585 ± 0.005, and M' 2 = 0.184 ± 0.006. The solid 
lines represent the estimated mean value for each quantity. 
These results are clear signs for a discontinuous transition 
[2411321 139 . as previously found for the random graph [341135] . 



the variables when the biggest jump occurs for each real- 
ization, and the maximum of M 2 . One observes that all 
three quantities are within error bars asymptotically in- 
dependent on the system size: J = 0.415 ± 0.005, Poo = 
0.585 ± 0.005, and M' 2 = 0.184 ± 0.006. These results 
suggest that the transition is discontinuous [Ml 1311 132] , 
as was previously reported for the random graph |341 135) . 
We remark that the numerical values of J, P^, and M' 2 
are within error bars identical to the ones found for a dif- 
ferent discontinuous percolation model in Refs. [Ml [52] . 

Unless indicated otherwise, on the square lattice, re- 
sults have been averaged over 10 8 samples for the small- 
est system size (32 2 ) and 70 samples for the largest one 
(8192 2 ) and error bars are smaller than the symbol size. 
Random numbers have been generated with the algo- 
rithm proposed in Ref. [10] • To keep track of the cluster 
properties as the fraction of occupied bonds increases, we 
have considered the labeling scheme proposed by New- 
man and Ziff [41j [42] , related to the Hoshen-Kopelman 
algorithm |4"5] . 

To determine the bond occupation fraction p c at which 
the percolation transition occurs, we measure for every 
sample the position p c j of the largest jump J, the posi- 
tion Pc.m of the maximum of M 2 , and the point p c _s at 
which spanning in one specific direction occurs. Extrap- 
olating these estimators to the thermodynamic limit, we 
find p c = 1.000 ± 0.002. In Fig. [5] we see p c — p c ^x as 
a function of the linear system size P, for all three esti- 
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FIG. 5. p c — Pc,x, with p c = 1.000 ± 0.002, as a function of 
the linear system size L for the three estimators, the position 
of the jump p C: j (o), the position p c ,M (*) of the maximum of 
M2, and the point p Ct s (□) at which spanning in one specific 
direction occurs. Asymptotically, we observe (p c — Pc.x) ~ 
L~ a with a = 0.5 ±0.2. 



mators. Asymptotically, the points follow a power law 
(Pc — Pc,x) ~ L~ a with a — 0.5 ± 0.2. In addition, with 
these results we can obtain a scaling plot for the standard 
deviation of the order parameter Xoo , 



Xc 
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where (•) is an average over independent realizations. A 
data collapse is obtained if Xoo(p, L) is plotted as a func- 
tion of (p — p c )L a , see Fig.[6j The fact that the maximum 
°f Xoo is independent on the system size is a further ev- 
idence for a discontinuous transition [2H|3U[5H|. Here, 
Xoo is defined as an intensive variable [see also Eq. 
as the standard deviation of the largest cluster size per 
lattice site. In the thermodynamic limit, this quantity 
is zero for continuous transitions and non-zero for dis- 
continuous transitions. For the considered system sizes, 
the peak height is independent of L. Results have been 
averaged over at least 1.5 x 10 3 samples. 

In Fig. [3] one observes that slightly below the percola- 
tion threshold the clusters of the BFW model are com- 
pact, but with a corrugated surface. This behavior has 
also been reported for another discontinuous percolation 
model, where the fractal dimension of the cluster perime- 
ter was found to be consistent with the one of the bridge 
line in bridge percolation and the watershed [24, 32, 44- 
[55] . Compact clusters with fractal surface have also been 
found in irreversible aggregation at high concentration 
[54] . For the BFW model, we measure at the stage where 
the number of clusters in the system equals two, for ev- 
ery sample the number A of bonds that would connect 
these two clusters if they would be occupied. The inter- 
face size [55] is observed to scale asymptotically with the 
system size, A ~ L d f (see Fig. [7]). To estimate the frac- 
tal dimension more accurately, we plot the local slopes 



FIG. 6. (Color online) Scaling plot for the standard deviation 
of the order parameter Xoo plotted as a function of the scaling 
variable (p — p c )L a , with p c = 1 and a — 0.5, for different 
linear system sizes L [32 (v), 64 (x), 128 (*), 256 (□), 512 (■), 
1024 (o), 2048 (•), 4096 (A), and 8192 (A)]. Note that this 
value of the exponent a agrees within error bars with the one 
obtained from the asymptotic behavior of (p c — p c ,x) ~ L~ a , 
a — 0.5 ± 0.2, see Fig.[5] One observes that Xoo is non-zero 
in the thermodynamic limit, as expected for a discontinuous 
transition. 



c?l = log( J 4i/ J 4 i / 2 )/l°g(2) as a function of the inverse 
system size L~ 5 , see lower right inset of Fig.[7j which 
yields an exponent of 1.49 ± 0.02. This result is consis- 
tent with the fractal dimension of the interface as ob- 
tained from box counting (see upper left inset of Fig.[7]) 
which yields an exponent of 1.49 ± 0.08. Interestingly, 
this value is different from the one of Ref. [24], where 
df = 1.23±0.03, but similar to the surface fractal dimen- 
sion in irreversible aggregation at concentration unity, 
where df = 1.60 ± 0.08 [54]. 

We note that for the BFW model clusters are com- 
pact and only the perimeter is fractal, with dimen- 
sion df. For 2D random percolation both the cluster 
mass and perimeter are fractals with dimension 91/48 = 
1.89583 ... p and df = 7/4 = 1.75 [55][56], respectively. 



We also simulated the BFW model on the 2D triangu- 
lar lattice and measured, in analogy to Fig. [4j the size of 
the jump J, the value of the order parameter imme- 
diately before the jump, and the maximum of M' 2 (not 
shown). Like on the square lattice, we found the three 
quantities to be asymptotically independent on the lat- 
tice size: J = 0.415 ± 0.005, P x = 0.585 ± 0.005, and 
M'2 = 0.183 ± 0.007, within error bars identical to the 
results for the square lattice. This is an evidence that 
the BFW percolation transition is discontinuous also on 
the triangular lattice. 
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FIG. 7. Size dependence of the number of interface bonds 
A (o). For large systems the interface size scales with the 
system size, A ~ L df . The straight line is a guide to the 
eye given by 0.824X 1 ' 47 . In the lower right inset we see the 
local slopes cLl = \og(AL/ A L /2)/ log(2) (•) as function of the 
inverse system size L -0 ' 5 . Extrapolating to L — > oo yields an 
exponent of 1.49 ± 0.02. Box counting results are shown in 
the upper left inset: we see the rescaled number of occupied 
boxes BL~ d f as a function of the linear box size n for two 
different system sizes: L = 8192 (■) and L = 4096 (□), the 
exponent is 1.49 ±0.08. Results have been averaged over 10 8 
samples for the smallest system size (16 2 ) and 10 3 samples 
for the largest one (8192 2 ). 




FIG. 8. (Color online) Fraction of sites in the largest cluster 
P x (solid red line), stage per site k/N (o), occupied bonds 
per step t/u (dashed upper blue line), and g(k) (dotted lower 
black line) as a function of the fraction of occupied bonds p, 
for the BFW model on the simple-cubic lattice. The system 
size is N = 64 3 , the measured values are averages over 10 3 
samples. In the inset we see the same variables, measured for 
a single realization. 




FIG. 9. (Color online) Fraction of sites in the largest cluster 
Poo as a function of the fraction of occupied bonds p on the 
simple-cubic lattice for different lattice sizes L (from left to 
right: 8, 16, 32, 64, 128, 256, and 512). 



III. BEHAVIOR IN THREE DIMENSIONS AND 
MULTIPLE JUMPS 



In Sec.[TT]we considered the BFW model on the square 
lattice and found evidence for a discontinuous transi- 
tion. We proceed to analyze the BFW model on the 
simple-cubic lattice. In Fig. [8] we see the analog of Fig. [I] 
for the simple-cubic lattice, namely, the fraction of sites 
in the largest cluster P^, stage per site k/N, occupied 
bonds per step t/u, and g(k) as a function of the frac- 
tion of occupied bonds p. One observes that the frac- 
tion of sites in the largest cluster P^ grows stepwise 
(see also Fig. [9]), which hints at the presence of multi- 
ple stable components, as reported for the random graph 
[331 [35] ■ As on the square lattice (see Fig.JT]), the size 
of the largest cluster is bounded by the stage k. The 
fraction of accepted bonds t/u decreases initially, ulti- 
mately reaching the asymptotic value of g(k) — 1/2 be- 
fore increasing slightly once only one cluster remains. In 
stages where s max is essentially constant, the fraction of 
accepted bonds increases due to the addition of redun- 
dant bonds. 

To analyze this effect in detail, we proceed as follows. 
For every sample we measure the Nj largest increases in 
the order parameter Poo, denoted as Jj, and the fraction 
of occupied bonds p c ,j.i at which they occur, and average 
over many samples. This definition contains the above 
ones of J and p c ,j, since J = J\ and p Ci j = p c ,j,i- Let us 
consider the six largest increases in the order parameter, 
i.e., Nj = 6. In the 2D square lattice, only one jump 
occurs: one observes that, extrapolated to L — > oo, the 
occupation fractions p c ,j,i are identical within their error 
bars, see Fig.[l0j[a), and the corresponding jump sizes Ji 
are macroscopic (not shown). This means that in the 
thermodynamic limit the fraction of bonds that is added 
between two macroscopic increases of the order parame- 
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FIG. 10. Bond occupation fractions p c ,j,i at which the 
six largest jumps in the order parameter occur for the BFW 
model as a function of the inverse system size, for (a) the 
square lattice and (b) the simple-cubic lattice. While for 
two dimensions all p c ,j,i can be extrapolated to the same 
value p c = 1.000 ± 0.002, in three dimensions four distinct 
jumps are visible within error bars for L — > oo (compare also 
Fig.|8|. For the square lattice, a = 0.5 ± 0.2 and the solid 
line is a guide to the eye. On the simple-cubic lattice we find 
p cA = 0.756±0.005, p c , 2 = 0.623±0.007, p c , 3 = 0.578±0.006, 
and p Ct i = 0.539 ± 0.008. For (b) results have been averaged 
over 10 8 samples for the smallest system size (8 3 ) and 40 sam- 
ples for the largest one (512 3 ). 



ter is zero. In this sense, for two dimensions, all jumps 
happen at the same instant p — p c . In contrast, as shown 
in Fig. 10 Id), in three dimensions, the macroscopic jumps 
occur at different bond occupation fractions p c ,j,i- This 
is consistent with multiple stable clusters emerging and 
coexisting when added bonds are mainly redundant ones. 
In Fig. [IT] we see the jump size J as a function of the frac- 
tion of occupied bonds p for the simple-cubic lattice. We 
divide the considered range of bond occupation fractions 
p into 512 bins of equal size. In each bin the maximum of 
the changes in the order parameter is measured and aver- 
aged over several samples. In the plot, one sees four peaks 
with heights increasing or being constant with increasing 
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FIG. 11. (Color online) Maximum change J in the order 
parameter in the BFW model on the simple-cubic lattice as 
a function of the fraction of occupied bonds p for different 
linear system sizes L [128 (□), 256 (o), and 512 (A): peaks 
become narrower with increasing L] . The arrows indicate the 
positions of the four largest jumps, extrapolated to L — > oo: 
p = 0.756, 0.623, 0.578, and 0.539 (see Fig. fTob . The range 



between zero and unity for p was divided into 512 bins of 
equal size and in every bin the maximum of the occurred 
changes in the order parameter was measured. Results have 
been averaged over 1.6 x 10 samples for the smallest system 
size (128 3 ) and 2.4 x 10 3 samples for the largest one (512 3 ). 



lattice size. Since these peaks are well-separated and do 
not seem to approach each other for increasing lattice 
sizes, this corresponds to four macroscopic jumps occur- 
ring at distinct values of p. This is consistent with the 
results shown in Fig. [To|b) and [8] the arrows in Fig. '. 
indicate the positions of the four largest increases in the 
largest cluster size, as obtained in Fig. 10 'b). We remark 



that the number of jumps measured at distinct fractions 
of occupied bonds depends on the resolution, which in 
Fig. [IT] is related to the bin size. An interesting possi- 
bility would be that the fraction of sites in the largest 
cluster Prx, as function of the fraction of occupied bonds 
p actually behaves like a Devil's Staircase, with an infi- 
nite hierarchy of jumps. In the next section, we study the 
tree-like version of the model where solely merging bonds 
are occupied, i.e., bonds connecting edges belonging to 
different clusters. 



IV. TREE-LIKE BFW MODEL 

The BFW model as defined in Sec. ITT1 seems to exhibit a 
discontinuous transition on the square and simple-cubic 
lattices. We found that in three dimensions multiple 
jumps of the order parameter occur at distinct bond oc- 
cupation fractions in the thermodynamic limit, which can 
only occur when redundant bonds are occupied. We now 
consider a modified version of the BFW model where re- 
dundant bonds are never occupied, i. e., bonds connecting 
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FIG. 12. (Color online) Fraction of sites in the largest cluster 
Poo as a function of the fraction of occupied bonds p for (a) 
the tree-like BFW model and (b) tree-like classical percolation 
on hypercubic lattices of dimension two to seven (from right 
to left with increasing dimension). For the BFW model the 
transition is abrupt. In the classical case, the transition is 
continuous and the thresholds are compatible with results for 
the number of clusters at the critical point 57 •">!) . The insets 
show the fraction of sites s/N in the three largest clusters [1st 
(□, upper curve), 2nd (A, middle curve), and 3rd (o, lower 
curve)] as a function of the fraction of occupied bonds p for 
single typical realizations of (a) the tree-like BFW model and 
(b) tree-like classical percolation on the hypercubic lattice of 
dimension seven of size 6 7 . While for the BFW model the 
largest cluster grows by overtaking, i.e., the second largest 
cluster overtakes the largest one which becomes the second 
largest one, for classical percolation the largest cluster grows 
directly by merging with smaller clusters. For the main plots, 
the system sizes L d are 2048 2 , 161 3 , 45 4 , 21 5 , 13 6 , and 9 7 for 
both models. Results have been averaged over at least 2 x 10 3 
samples. 



sites that are part of the same cluster arc not sampled 
and, therefore, all clusters are tree-like. This version was 
also considered for random graphs in Ref. 34; under the 
name restricted BFW model. 

Figure [12] shows as a function of the fraction of 
occupied bonds p for the tree-like BFW model and, for 



the sake of comparison, we also include the results for 
tree-like classical percolation. We observe that under 
this constraint, for the BFW model the transition seems 
abrupt in all considered dimensions, whereas in the clas- 
sical case the transition is continuous and the percolation 
thresholds are consistent with the values expected for the 
number of clusters at the critical point (see Ref. [571155] 
and references therein). In particular, at p c , the num- 
ber of clusters per site (N c /N) in tree-like percolation 
is N c /N = 1 — p c d, which on the square lattice, yields 
p c w (7 - 3^3) /4 [57H59]. The insets in Fig.[l| show 
the fraction of sites s/N in the three largest clusters as 
a function of the fraction of occupied bonds p for single 
realizations of the tree-like BFW model [Fig.[l2^a)] and 
tree-like classical percolation [Fig.|12[b)]. One observes 
that for the BFW model the size of the largest cluster 
grows due to two smaller clusters merging together to 
become the new largest (growth by overtaking), while 
in the classical case the largest component mainly grows 
by adding small pieces [35]. In the inset of Fig.[l2|a) 
the symbols corresponding to the largest cluster form the 
curve on top. When a jump occurs, the size of the largest 
cluster changes by overtaking: for example at the second 
last jump shown in the inset, the second largest clus- 
ter merges with the third largest cluster and becomes 
the largest cluster, while the previous largest cluster is 
then the second largest one. For the complete graph it 
has been shown that all significant growth of the largest 
cluster size is due to overtaking [35]. 

We now analyze the behavior of the tree-like BFW 
model on hypercubic lattices. To estimate the percola- 
tion thresholds for the infinite system, we consider the 
estimator pj, i. e., the bond occupation fraction at which 
the largest change in the order parameter occurs (see 
Sec.[n]) and study its size dependence. For dimensions 
two to seven we obtain, respectively, p c = 0.500 ± 0.001, 
0.333±0.002, 0.248±0.002, 0.198±0.003, 0.165±0.004, 
and 0.141 ± 0.009 (see Fig. 13). These values are within 
error bars equal to 1/d. 



In Ref. [32] it was shown that 
for any percolation model on hypercubic lattices of di- 
mension d with finite number of clusters at the thresh- 
old p c > 1/d. If, in addition, the models are tree-like, 
p c = 1/d was shown. The values of the jump J are 
summarized in Fig. |14| where J is plotted as a function 
of the system size N for the tree-like BFW model in all 
considered dimensions. One observes that the jump value 
converges towards a constant value for large system sizes, 
which is a strong evidence of a discontinuous transition 

EH]. 

In the following we investigate the cluster-size distri- 
bution at the threshold for high dimensions. In Fig. 15 
we see a plot of the number of clusters of size s per site n s 
as a function of the cluster size s for hypercubic lattices 
of dimension three, five, and seven. To reduce statistical 
fluctuations, we considered the accumulated number of 
clusters N s , defined as N s — Y^r=s^ n r- K can be ob- 
served that the cluster-size distribution is composed of 
two parts as the dimension increases, in agreement with 
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FIG. 13. (Color online) Percolation threshold estimator p c — 
p Ct j as a function of the inverse system size N~ for the tree- 
like BFW model on the hypercubic lattice of dimensions two 
to seven. Extrapolating to L — > oo, we findp c = 0.500±0.001, 
0.333 ±0.002, 0.248 ±0.002, 0.198 ±0.003, 0.165 ±0.004, and 
0.141 ± 0.009. Note that p c seems to be 1/d. In the inset we 
see the same data for p c —Pc,j, plotted as a function oiN = L d 
in double-logarithmic scale, where the solid line is a guide to 
the eye of the form 0.018-/V -0 ' 7 . Results have been averaged 
over at least 7.8 x 10 6 samples for the smallest system size 
(64 2 , 16 3 , 8 4 , 5 5 , 4 6 , and 4 7 ) and over at least 100 samples 
for the largest one (8192 2 , 512 3 , 91 4 , 37 5 , 20 6 , and 16 7 ). 
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FIG. 15. (Color online) Accumulated cluster-size distribu- 
tion N 3 as a function of the cluster size s for the tree-like BFW 
model on hypercubic lattices of dimension three (■,□), five 
(A, A), and seven (T,V). N s = J^fT - , 1 n r , where n s is the 
number of clusters of size s per lattice site N — L d , exclud- 
ing the largest cluster of size s max . One observes that the 
distribution falls off sharply for small cluster sizes and has a 
region on the right, qualitatively similar to the behavior on 
the random graph 34, 35]. For comparison, the inset shows 
N s as a function of s for tree-like classical percolation on the 
square lattice at p c = (7-3^3)/4 57-59!. Results are shown 
for system sizes of 8192 2 (o) and 16384 2 (•). Asymptotically, 
the points follow a power law, N s ~ s 1_T , where r = 187/91. 
All results have been averaged over 100 samples. 
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FIG. 14. (Color online) The largest change in the order 
parameter: J is plotted as a function of N — L d for the tree- 
like BFW model on hypercubic lattices of dimension two to 
seven. Asymptotically, the jump is independent on system 
size, indicating a strongly discontinuous transition. The solid 
lines are guides to the eye indicating the following constant 
values: 0.415±0.005, 0.414±0.002, 0.412±0.003, 0.410±0.003, 
0.409 ±0.002, and 0.409 ±0.002, for dimension two to seven. 



the behavior on the random graph 34, 35 . The number 
of small clusters at the threshold increases with dimen- 
sion. Thus the cluster-size distribution for the tree-like 
BFW model behaves in a different way than in classical 
percolation at the critical point, where asymptotically a 



power law holds, 
Fig.pl). 



or N x 



(see inset of 



Besides finding some evidence that the tree-like BFW 
model on the hypercubic lattice exhibits a discontinu- 
ous transition, the results presented in this section also 
show that plateaus in the evolution of the fraction of 
sites in the largest cluster on the simple-cubic lat- 
tice (see Fig. [8]) are due to redundant bonds, which are 
not present in the tree-like case, in agreement with what 
has been reported for the random graph [341 135] . How- 
ever, on the two-dimensional square lattice, the jumps 
occur at one unique fraction of occupied bonds, irrespec- 
tive of whether the tree-like version is considered or not. 
This result is in agreement with the results of Moreira 
et ai, showing that to obtain a discontinuous percola- 
tion transition in two dimensions, the control of inter- 
nal/redundant bonds is irrelevant |30j . 
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V. FINAL REMARKS 

We analyzed the BFW model [34H36] on the lattice 
and found numerical evidence for a strongly discontin- 
uous transition 38J in two and three dimensions. The 
surface fractal dimension of the BFW model in two di- 
mensions was measured to be dt = 1.49 ± 0.02. Interest- 
ingly, this is different from the fractal dimension reported 
for another discontinuous percolation model, 1.23 ± 0.03 
[241 153| , but similar to the one found in irreversible aggre- 
gation at high concentration, 1.60 ±0.08 [54] ■ The BFW 
model on the simple-cubic lattice was found to exhibit 
multiple distinct jumps of the largest cluster size in the 
thermodynamic limit. We further analyzed the tree-like 
BFW model on hypercubic lattices of dimension two to 
seven, finding that the jump size is asymptotically inde- 
pendent on system size and determined the percolation 
thresholds. These results imply that the coexistence of 
multiple stable components in the BFW model is due to 
redundant bonds. This can be seen from the fact that 
the tree-like version of the BFW model only exhibits 
one jump while multiple jumps at distinct fractions of 
occupied bonds appear only for the simple-cubic lattice 
when redundant bonds are also sampled. We also consid- 
ered the evolution of the cluster-size distribution at the 
threshold with increasing dimension. 

The results presented here are consistent with the be- 
havior of the BFW model on the random graph [M] [35] . 
The mechanisms leading to discontinuous percolation in 
this case are compatible with the required ones according 



to Ref. [30] . namely, the homogenization of cluster sizes 
and the control of merging bonds. Identifying percola- 
tion models with discontinuous transitions is interesting 
since the original proposal [3J and related models have 
been shown, both analytically and numerically, to yield 
continuous transitions in the thermodynamic limit [19j . 
Nevertheless, these initial models are characterized by 
a macroscopic jump for finite networks even with 10 18 
nodes [201111]. Note that, the Internet and the WWW 
have less than 10 12 nodes. 

It would be interesting to understand the finite-size 
scaling of discontinuous percolation models in a more 
rigorous analytical manner since, in the presence of re- 
dundant bonds, it differs from the one of equilibrium 
first-order phase transitions [39]. In addition, interesting 
fractal dimensions emerged in the investigation of cluster 
perimeters in discontinuous percolation [3H |4"4"H5"2"] which 
could also be studied for the BFW model in dimensions 
higher than two. 
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